Limits on the Curie temperature of (III,Mn)V ferromagnetic semiconductors 
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Mean-field-theory predicts that the Curie temperature T c 
of a (III,Mn)V ferromagnet will be proportional to the va- 
lence band density-of-states of its host (III,V) semiconductor, 
suggesting a route toward room-temperature ferromagnet ism 
in this materials class. In this Letter, we use theoretical es- 
timates of spin-wave energies and Monte-Carlo simulations 
to demonstrate that long-wavelength collective fluctuations, 
neglected by mean-field theory, will limit the critical temper- 
ature in large density-of-states materials. We discuss impli- 
cations for high T c searches. 

The recent discovery Q of ferromagnetism at relatively 
high temperatures (T c > 100K) in III-V compound semi- 
conductors containing Mn has generated intense inter- 
est, mainly because of the new technological roadways 
that would be opened by room temperature ferromag- 
netism in semiconductors with favorable materials prop- 
erties. The search for systems in this materials class with 
higher critical temperatures is an important current ac- 
tivity that has been guided thus far by mean-field the- 
oretical l^ll considerations. In this Letter we address 
the importance of collective magnetization fluctuations, 
neglected by mean-field theory, in limiting the critical- 
temperature and discuss the implications of these con- 
siderations for high-T c searches. 

Our analysis is based on the kinetic-exchange model for 
interactions between the Mn spins and band electrons, 



H — Hq + J, 
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where Si describes a Mn spin with spin length S = 5/2 at 
site Ri, s(r) is the band-carrier density, and J p d > rep- 
resents the exchange integral. Hq represents a simplified 
single-parabolic-band model for the host semiconductor 
valence bands. 

The simplest treatment of this model is a mean field 
theory which takes the magnetizations of carriers and ion 
spins to be uniform in space and neglects correlations 
between them. A straightforward calculation yields for 
the critical temperature in mean field approximation 
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where g* is the g-factor of the carriers and xp 1S their 
Pauli susceptibility, which is proportional to the effective 
band mass. This observation has given rise to concrete 



predictions for critical temperatures for several host semi- 
conductors based on their different band masses [g| . 

In mean-field theory, ferromagnetism occurs because 
the penalty in entropic free energy paid to polarize the 
Mn spins vanishes at T = 0. Any coupling to a band- 
electron system with a finite magnetic susceptibility is 
sufficient to yield ferromagnetism. While this mean- 
field theory probably captures much of the physics of 
(III,Mn)V ferromagnetism, it has a qualitative deficiency 
which will have an important quantitative impact on T c 
predictions in circumstances we identify below. Mean- 
field theory fails to account for the small energy cost of 
magnetization configurations in which spin orientations 
vary slowly Q in space, reducing the average magnetiza- 
tion but maintaining local correlations between Mn and 
band-electrons spin orientations. In the following para- 
graphs we estimate the critical temperature for the case 
when these collective excitations dominate thermal mag- 
netization suppression. 

Isotropic ferromagnets Q have spin-wave Goldstone 
collective modes whose energies vanish at long wave- 
lengths, 



E{k) = Dk 2 
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where k is the wavevector of the mode. Each spin-wave 
excitation reduces the total spin of the ferromagnetic 
state by 1. The coefficient D is inversely proportional 
to the saturation magnetization and proportional to the 
exchange constant A of classical micromagnetic theory, 
that parameterizes the free-energy cost of spatial varia- 
tions in magnetization orientation. We have previously 
presented a theory of spin-wave excitations in (III,Mn)V 
ferromagnets Q . These collective excitations are not ac- 
counted for in the mean-field approximation. If the spin 
stiffness is small, they will dominate the suppression of 
the magnetization at all finite temperatures and limit the 
critical temperature. A rough upper bound on the re- 
sulting critical temperature, T£ oU , can be found by using 
the T = stiffness value and finding the temperature 
where the number of excited spin waves per unit volume 
equals the spin per volume of the ground state, 
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A Debye wavevector 0] k D = (fo^N) 1 ^ cuts off the 
sum over wavevectors at the correct number of modes, 
and n{E) is the Bose occupation number. We find that 
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for 5 > 5/2. To obtain this equation, we have as- 
sumed that the spin waves can be approximated as non- 
interacting Bose particles, replaced the dispersion by the 
long- wavelength limit Eq. (||) , and noted that the critical 
temperature estimate is proportional to Dk 2 Dl justifying 
the use of the classical expression for the mode occupa- 
tion number n(k) « kBT/E(k) — 1/2. These consider- 
ations set an upper bound on the critical temperature 
which is proportional to the spin-stiffness, a bound not 
respected by the mean-field theory. A familar example 
of ferromagents in which long range order is suppressed 
by long-wavelength collective excitations is provided by 
the ferromagnetic transition metals Fe, Co, and Ni. In 
that case, an expression similar to Eq. (^|) and propor- 
tional to the micromagnetic exchange constant, predicts 
critical temperatures with 20% accuracy |^|| whereas 
mean- field-theory overestimates T c by factors of 5 to 10. 
In (III,Mn)V ferromagnets, we will see that both mean- 
field and collective regimes can occur depending on car- 
rier density and host semiconductor band parameters. 
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FIG. 1. Critical-temperature-limit regimes. In the 
mean-field regime T c is limited by individual Mn spin fluc- 
tuations. In the collective regimes, the critical temperature is 
limited by long-wavelength fluctuations with a stiffness pro- 
portional to the bandwidth for weak (RKKY) exchange cou- 
pling and inversely proportional to the bandwidth for strong 
exchange coupling. At the solid line T C MF = T c c °". Dashed 
lines: expansions for large and small A/cf, Eqs. (^) and (Q), 
and the crossover from the RKKY to the strong coupling col- 
lective regime. 

Our theoretical results for (III,Mn)V ferromagnet col- 
lective modes lead to physically transparent results for 
the spin-stiffness in both strong and weak exchange cou- 
pling limits. The dimensionless parameter which char- 
acterizes the strength of the exchange coupling is the 
ratio A/ep, where ep is the band-system Fermi en- 
ergy and A = J p dNS is its mean-field spin splitting 



at T = 0. For small A/e F , the RKKY regime f|, ex- 
change coupling is a weak perturbation on the band sys- 
tem. In this regime we find that D = 6/(12k F ) where 
S = Jpd</2 = (3/8S)(n/N)(A 2 /e F ) is the energy cost 
of an uncorrelated spin reduction at a single Mn site. 
Note that in this regime 6 ~ T^ F and that 
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In the weak coupling regime mean-field theory is reliable 
only for n/N <C 1, as expected since in this case the 
RKKY interaction has a range which is long compared 
to the distance between Mn spins. In the large A/ep 
regime exchange coupling completely polarizes the band- 
electron system. In this case (and for n <C 2NS) we find 
that D = (n/2NS)(e F /k F ). For a fully polarized band 
the energetic cost of varying the moment orientation di- 
rection is entirely due to band kinetic energy. We, thus, 
obtain as a third T r bound 
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The different regimes deduced from these considerations 
are illustrated in Fig. |l|. 

To substantiate these qualitative considerations, we 
have performed Hybrid-Monte-Carlo pj simulations, 
treating the Mn spins as discrete classical degrees of free- 
dom, an approximation that is justified near the critical 
temperature. We allow for disorder by choosing the Mn 
positions randomly. Microscopic p — d exchange physics 
is modeled by allowing the interaction to have a finite 
range fill . We simulate this by replacing the delta 
function by a Gaussian distribution in Eq. ([!]). 

An exhaustive description of our Monte Carlo ap- 
proach, including a detailed account of all technical as- 
pects such as thermalization procedures, finite-size effctes 
etc., will be given elsewhere || . Here we shall, for brevity, 
concentrate on the results. In the following we consider 
the strong coupling regime where mean-field theory is 
not reliable and finite-size effects in our simulations are 
small. One important finding is that randomness in the 
Mn positions can enhance the spin stiffness (for large J p d 
by up to factor of two) in comparison to the estimates 
based on course-grained Mn spins with Debye cutoffs 
discussed above. The importance of disorder was also 
emphasized by Wan and Bhatt [Q for a model of in- 
teracting magnetic polarons where both the ion and the 
carrier spins are treated classically, while we account for 
the quantum- mechanical nature of the free carriers in our 
simulation. 

Figure ^ shows a typical magnetization curve for man- 
ganese ions and carriers magnetizations in the collective 
regime. The maximum of the finite-size Mn susceptibil- 
ity marks the ferromagnetic transition temperature. The 
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lower inset compares the global and local free carrier po- 
larizations, mioc = (\s(r)\/n(r)} . As anticipated by our 
qualitative discussion, the typical free carrier local band 
polarization remains large above the critical temperature; 
it would vanish in a mean-field picture. Ferromagnetism, 
and the technologically useful robust collective physics it 
gives rise to, disappears in this case only because of the 
loss of long-range spatial coherence. 
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FIG. 2. Magnetization curves for Mn ions and carriers. 
The upper inset shows the magnetic susceptibility due to the 
ions, while in the lower inset the local and the global spin 
polarization per carrier is plotted. The data was obtained for 
a system of 540 Mn spins and 54 carriers. 
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FIG. 3. The critical temperature T c as a function of the 
carrier mass (left panel) and as a function of the carrier den- 
sity (right panel). The results of the Monte Carlo runs are 
compared with the mean-field predictions. 

The left panel of Fig. || shows the critical temperature 
as a function of the carrier effective mass. Mean-field 
theory predicts that T c will grow linearly with increasing 
mass. The Monte-Carlo results, however, are substan- 
tially lower and show a saturation of T c at carrier masses 



close to the bare electron mass. For even higher masses, 
we expect T c to decrease, reflecting the reduction of spin 
stiffness expected in this regime. 

The right panel of Fig. [}] compares the mean-field pre- 
diction and the Monte Carlo results for T c as a function of 
the carrier density. For higher carrier densities we expect 
ferromagnetism to give way to spin-glass order. 

We conclude from the present work that high critical 
temperatures cannot be achieved simply by narrowing 
the free carrier band or placing its Fermi energy at a 
density-of-states peak in order to enhance its Pauli mag- 
netic susceptibility \p- It will a ls° be necessary to en- 
gineer a suppression of collective magnetization fluctua- 
tions. 

We thank Glenn Martyna for pointing out to us the 
striking advantages of the Hybrid Monte Carlo algorithm, 
and T. Dietl, T. Jungwirth and H. Ohno for useful dis- 
cussions. We acknowledge support from the Deutsche 
Forschungsgemeinschaft, the National Science Founda- 
tion, and the Indiana 21st Century Fund. 



[1] 



[2] 



For reviews see H. Ohno, J. of Magnetism and Magnetic 
Materials 200, 110 (1999), H. Ohno, Science 281, 951 
(1998); for a more general introduction to the emerging 
field of spintronics see G. A. Prinz, Science 282, 1660 
(1998). 

T. Dietl, A. Haury, and Y. M. d'Aubigne, Phys. Rev. 
B 55, R3347 (1997); T. Dietl et ai, Science 287, 1019 



(2000); T. D ietl, H. Ohno, and F. Matsukura, cond 



mat/0007190 



[3] 



[4] 



[5] 



[6] 



[7] 



T. Jungwirth, W. A. Atkinson, B. H. Lee, A. H. MacDon- 
ald, Phys. Rev. B 59, 9818 (1999); B.H. Lee, T. Jung- 
wirth, A.H. MacDonald, Phys. Rev. B 61, 15606 (2000); 
M. Abolfath, T. Jungwirth, J. Brum, A. H. MacDonald, 
Phys. Rev. B 63, 054418 (2001). 

For these qualitative considerations we replace the dis- 
crete Mn spins by a continuous distribution. This is a 
good approximation when the Mn density is much larger 
than the free-carrier density. 

Because of spin-orbit coupling in the valence band the 
magnetic anisotropy energy is non-zero in (III,Mn)V fer- 
romagnets and there is a resulting gap in the spin-wave 
excitation spectrum. [J. Konig, T. Jungwirth, H.H. Lin, 
A.H. MacDonald, to be submitted to Physical Review B]. 
However, this gap is quite small and does not change the 
essence of the present discussion. 

J. Konig, H.H. Lin, a nd A.H. MacDonal d, Phys. Rev. 



Lett. 84, 5628 (2000); |cond-mat/001047l| , to appear in 
"Interacting Electrons in Nanostructures" , edited by R. 
Haug and H. Schoeller (Springer 2001). 
At finite temperature interactions between spin waves 
will reduce the spin stiffness. To model this effect a self- 
consistent spin-wave scheme would be required. Here, 
however, we use the zero-temperature spin stiffness for 



3 



an upper bound. 
[8] J. Schliem ann, J. 



mat/0012233 



Konig, A. H. MacDonald, cond 



[9] O. Gunnarsson, J. Phys. F 6, 587 (1976). 
[10] S. Duane, A. D. Kennedy, B. J. Pendleton, D. Roweth, 

Phys. Lett. B 195, 216 (1987). 
[11] A. K. Bhattacharjee, C. Benoit a la Guillaume, Solid 

State Comm. 113, 17 (2000). In the present Monte Carlo 

simulations we use ao — O.lnm. 

[12] X. Wan, R. N. Bhatt, |cond-mat/0009161 . 



4 



